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1.  Controllability  of  infinite  dimensional  bilinear  systems 

A  standard  equation  governing  small  transverse  deflections  of  a  beam  is 

utt  +  uxxxx  +  P(t)v»xx  =  0  0  <  x  <  1  (1.1) 

where  u(x,t)  denotes  the  deflection,  p(t)  represents  the  axial  load  (a  real 
valued  function).  A  natural  control  problem  one  might  consider  is  the 
following:  Given  initial  data 

u(x,0)  =  u0{x)  ,  ut(x,0)  ®  (x)  (1.2) 

consistent  with  some  reasonable  set  of  boundary  conditions  and  terminal  data 

u(x,T)  *  uT(x)  ,  ut(x,T)  =  u'T(x)  (1.3) 

find  a  control  p(t)  e  L2(0,T]  which  accomplishes  this  transfer. 

For  example,  if  we  consider  the  cased  of  hinged  boundaries  u  *  uxx  =  0 
at  x  =  0,  1  and  expand 

u(x,t)  a  l  un(t)$n(x)  (1.4) 

u=1 

where  ♦n(x)  «  /2  sin  nux.  We  see  the  {un(t)}  satisfies  the  infinite  system 
of  ordinary  differential  equations 

un(t)  +  n4*4un(t)  -  p(t)n2w2un(t)  -  0  ,  (1.5) 

n  «  1,2,3,...  .  Mow  our  goal  of  state  transfer  can  be  viewed  as  desiring  one 
control  p(t)  which  transfers  each  Fourier  coefficient  of  the  initial  data  to 
the  desired  Fourier  coefficient  of  the  terminal  data.  This  problem  is  far 
from  trivial  since  (i)  we  are  attempting  to  control  an  infinite  number  of 
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modes  with  one  control  p(t),  (ii)  the  input-output  relation  p  ♦  {un}  is 
nonlinear. 

For  problems  of  the  type  where  the  inodes  decouple  as  in  (1.5),  Ball, 
Marsden,  and  Slemrod  [1]  developed  a  theory  which  identified  attainable 
terminal  data  for  a  reasonable  class  of  controls.  In  research  done  this  past 
year,  Slemrod  [2]  has  extended  these  results  to  a  class  of  control  systems  for 
which  the  inodes  do  not  decouple.  For  example,  Slemrod  theory  allows 
identification  of  attainable  terminal  data  for  control  systems  of  the  form 


un ( t ) 4  Xnun (t)  +  p(t)  J  Bjnn Xmum ( t )  *  ^ 

m=1 


(1.6) 


where  the  coupling  coefficients  B^  satisfy  special  hypotheses.  Slemrod 
applied  his  results  to  a  class  of  wave  equations  with  mixed  boundary 
conditions.  The  main  technique  was  to  perform  a  change  of  variables  that 
allowed  application  of  the  classical  Fredholm-Riesz-Schauder  theory  for  a 
linearized  problem  near  p  =  0  and  then  application  of  "local  onto"  theorem 
for  obtaining  the  nonlinear  result.. 

Future  research  will  be  based  on  this  approach,  possibly  employing  the 
Leray-Schauder  degree  theory. 

2.  Singular  perturbations  approach  to  optimization  problems. 

In  this  research  Slemrod  (in  collaboration  with  J.  Carr  of  Heriot-Watt 
University,  Edinburgh,  U.K.  and  M.  E.  Gurtin,  Carnegie-Mellon  Univ. , 
Pittsburgh,  Pa.),  we  are  considering  the  following  problem  motivated  by  the 
equilibrium  theory  of  phase  transitions. 


T 


Find  a  function  w(x),  u  <  x  <  1,  so  that  w(x)  minimizes  the  "cost"  or 
"energy"  functional 

1 

E(w)  =  /  W(w'(x)dx 
0 


with  prescribed  boundary  conditions.  If  W  is  convex,  a  minimum  always  exists. 
However,  if  W  is  non-convex,  typically  only  a  "relaxed"  solution  exists,  i.e., 
a  solution  to  the  problem  exists  for  W  replaced  by  its  convex  hull.  In  our 
research  we  have  studied  the  singularly  perturbed  problem  of  minimizing 


1 

/  ew'*2  +  W(w'(x)dx  ,  e  >  0 

0  e  e 


subject  to  boundary  conditions.  We  have  completely  analyzed  the  effect  of 
such  higher  gradients  on  the  solutions  for  small  e  >  0.  In  fact,  we  have 
shown  that  introduction  of  the  e  term  appreciably  changes  the  situation  in 
that  solutions  to  the  second  problem  are  essentially  unique  and  have  unique 
limits  as  £  ♦  0.  In  this  sense,  incorporating  the  physically  expected  e  term 
in  the  cost  and  then  passing  to  the  limit  produces  an  essentially  unique 
solution  whereas  studying  the  relaxed  problem  immediately  produces  an  infinity 
of  solutions.  We  view  our  approach  as  a  natural  "admissibility"  criteria  for 
non-convex  cost  problems. 

One  paper  has  already  been  completed  [3].  A  second  is  in  preparation. 

3.  Nonlinear  continuum  mechanics  and  phase  transitions. 


In  the  study  of  phase  transitions  in  fluids  a  typical  model  given  in 
classical  thermodynamics  is  provided  by  the  equilibrium  configurations  of  a 
van  der  Waals  fluid.  Such  equilibria  show  the  fluid  may  exist  in  two  phases. 
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liquid  and  vapor.  Relatively  untouched  in  classical  thermodynamics  is  the 
non-equilibrium  case,  i.e.,  how  the  dynamics  of  phase  transitions  proceed. 
The  goal  in  this  research  was  to  study  precisely  this  subject  from  continuum 
mechanics  (not  statistical  mechanics)  point  of  view. 

While  earlier  work  has  dealt  with  the  "structure"  of  interphase  shocks, 
Slemrod  has  lately  been  considering  the  role  of  viscosity  and  capillarity 
effects  on  the  full  hydrodynamic  equations  of  mass,  momentum,  entropy,  and 
energy.  Specifically,  he  has  been  considering  the  problems  of  (i)  whether 
inviscid  hydrodynamics  is  really  a  good  model  for  nonequilibrium  phase 
transitions  and  (ii)  how  should  we  compute  numerically  solutions  to  the 
hydrodynamic  equations .  Tentative  answers  to  these  problems  seem  to  be 
(i)  inviscid  hydrodynamics  is  not  a  good  model  for  nonequilibrium  phase 
transitions  since  omitted  viscous -capillarity  terms  are  crucial  to  the 
limiting  results  and  (ii)  a  naive  Lax-Friedrichs  scheme  may  be  a  reasonable 
way  of  performing  isothermal  computations. 

In  related  work  with  J .  K.  Hunter  [4],  Slemrod  has  considered  dynamic 
phase  transitions  in  polymers.  Specifically,  they  have  proposed  a 
phenomo logical  explanation  for  the  appearance  of  the  "ripple"  phenomena  in 
melt  fracture  of  certain  polymers.  Their  main  idea  was  to  consider  a 
nonlinear  monotone  rate  of  strain  history  functional  for  the  shearing  stress 
in  an  isotropic  viscoelastic  fluid.  This  led  them  to  consider  the  initial 
history  valued  problem 

OB 

vt(x,t)  =  3  /  e”*3  o(vx(x,t-s ) )ds  +  y  (3.1) 

lx  0 

v(x,t)  -  vg(x, t)  ,  -•  <  t  <  0  , 


t 
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where  v(x,t)  is  the  fluid  velocity  at  a  point  x  at  time  t,  o  is  a  non-monotone 
function,  y  is  the  applied  pressure  gradient.  Slemrod  and  Hunter  showed  that 
solutions  of  (3.1)  exhibit  hysteresis  and  phase  separation  in  a  manner  similar 
to  observed  experimental  results. 

4.  Optimal  control  problems  arising  in  robotics. 

Slemrod  supervised  a  Master  in  Applied  Mathematics  project  of  Gabriel 
Amazigo  on  modeling  and  control  of  a  deformable  arm.  In  this  project  Amazigo 
first  used  standard  Timoshenko  beam  theory  to  obtain  a  set  of  ordinary 
differential  equations  governing  the  lowest  modes  of  an  axially  and  flexurally 
deformable  arm.  In  the  second  part  of  the  project  Amazigo  used  the 
Pontrayagin  maximum  principle  to  find  an  optimal  control  to  bring  the  arm  from 
a  deformed  state  back  to  its  state  of  rest.  This  work  is  currently  being 
written  up  by  Amazigo. 
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